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ABSTRACT 

We present, using a novel technique, a study of the angular distribution of satellite 
galaxies around a sample of isolated, blue host galaxies selected from the sixth data 
release of the Sloan Digital Sky Survey. As a complement to previous studies we 
subdivide the sample of galaxies into bins of differing inclination and use the systematic 
differences that would exist between the different bins as the basis for our approach. 
We parameterize the cumulative distribution function of satellite galaxies and apply 
a maximum likelihood, Monte-Carlo technique to determine allowable distributions, 
which we show as an exclusion plot. We find that the allowed distributions of the 
satellites of spiral hosts are very nearly isotropic. We outline our formalism and our 
analysis and discuss how this technique may be refined for future studies and future 
surveys. 
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1 INTRODUCTION 

The standard cosmological model that assumes cold dark 
matter particles and a cosmological constant as domi- 
nant components of our universe (LCDM) has reached a 
stage where cosmological tests at the scale of galaxies are 
possible (e.g. 



. .Ayila- Reese et alj|1998t van den Bosch 19981: 
Courteau fc Rbd Il99£t iGnedin et all booet iPizagno et all 
20071 ). Although powerful, these types of cosmological con- 



straints are hampered by uncertainties in our understanding 
of baryonic evo lution. This limitation has led to contr over- 
sial results fe.g. lMoore et al.lll999bl : lKlypin et al.ll999al ) and 
demonstrate the need for modifications to the model and/or 
for more detailed comparisons between theory and obs erva- 
tions (e.g. IValenzuela et al.ll2007l ; ISimon fc Gehall2007l ). 

Under the LCDM picture of structure formation, galac- 
tic sized halos are assembled by accret i ng smaller structures 
or subhalos rtKauffmann et al. | Il993l : iGhigna et al. I 1 19981 : 
iKlypin et al.lll999bl ). Consequently, the presence of satellite 
galaxies is a natural prediction of hierarchical models and 
provide important constraints on galaxy-scale dark matter 
halos that are not always accessible by other means. For 
example, it is well known that the radial distribution of 
satellites can constrain the radial struc ture of the larger, 
host dark matter halo (|Chen et al.ll2006l ). The satellite an- 
gular distribution may also provide valuable information be- 
cause there are indications that the subhalo population is a 
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good tracer of dark matter halo shape (|Zentner et al.ll2005l : 
lAgustsson fc Brainerdll2006l ). 

The pionering observational stu dy on the angular distri- 
bution of satellites was presented bv lHolmber 1 ljl974 ). This 
study found an excess of satellites along a direction perpen- 
dicular to the host galaxy disk . The result was later con- 
firmed by (jZaritskv et al]|l997l ). Since the generic LCDM 
halo is triaxial, a combination of orbital precession and 
anisotropic dynamical friction mig ht create such a "Holm- 
berg" effect |Penarrubia et akll2004l ). A considerable amount 
of subsequent work discusses the statistical significance and 
the observational syste matics of the result as well as differ- 
ent interpretations (e.g. iBailin et al]|2007l ). 

Another prediction of hierarchical models that could 
contribute to an anisotropic distribution of satellite galax- 
ies is the anisotropic accretion of material along large scal e 
filaments dPenarrubia fc Bensorj|2005l : IZentner et alil2005l ). 
Both the halo shape and its accretion history can imprint 
a signature on the satellite distribution with respect to the 
host galaxy and on the orientation of both the host and 
its satellites with respect to the surrounding large scale 
structure (e.g Aragon-Calvo et all 120071 ; lHahn etal 1 120071 ; 
IWang et al.1l2007l) . 

Several recent studies aimed to increase the statistical 
significance of the measurements and to discuss the system- 
atic effects with the hope of constraining t h e shap e of galaxy- 
scale dark matter halos. IKroupa et al.l ([2005 ); Ha rtwick 



(2000) discuss the case of the Mil ky Way. iKoch fc Grebel 
(2006) discuss the M31 system, and lMetz et all (|2007l ) diss- 



© 0000 RAS 



2 Steffen and Valenzuela 



cuss both. In these cases an anisotropic distribution was 
found, with relatively high significance, that favors a po- 
lar alignment. Simultaneously, a number of studies adopted 
the approach of selecting galaxies and their satellites from 
large surveys and obtain the opposite result; satellites pref- 
erentially align with t h e observed lon g axis of th e host 
llSales fc Lambasl 



I200i iBrainerdl [20051; 



Yang et alj 120061: 



Azzaro et alj 120071 ; iFaltenbacher et all 120071 ; Kang et al.l 
200^ 



TBailin et all (|2007T ) recently reviewed the situation and 
suggested that systematics related to the inclusion of non- 
relaxed systems (groups or clusters) precluded the detection 
of the Holmberg effect for late type galaxy hosts in those 
studies. The existing discrepancies indicate a need for larger 
samples, a better understanding and control of systematic 
effects, and new or complementary analysis approaches in 
order to have a clear understanding and robust interpreta- 
tion of the observations and to properly connect the obser- 
vations to theoretical models. 

1.1 Motivation 

Here we outline a new, complimentary method to study the 
distribution of satellite galaxies that may serve as a link 
between numerical work and observations and which may 
help to resolve some of the outstanding issues in this field. 
One challenge that exists in relating numerical studies to 
observations is centred on the system of coordinates that is 
used. From numerical simulations it is relatively straight- 
forward to identify the important physical axes of a galaxy 
scale dark matter halo. However, the relative orientation of 
the dark matter halo to the baryons is not always clear. 
In addition, observational studies of satellite galaxies use 
coordinate systems based upon the profile of the luminous 
matter as opposed to a system based upon the unobserved 
dark matter. 

This disconnect is particularly relevant if the distribu- 
tion of satellite galaxies is used as a tracer of the overall 
shape of the dark matter halo. Since observations only show 
a two-dimensional projection of the material, it is not possi- 
ble to identify which axes of a triaxial halo you are observing 
and how the system is inclined with respect to the line of 
sight. A particular ellipsoid can look either circular or highly 
elliptical depending upon the location of the observer. Thus, 
the observationally motivated coordinate system based upon 
the distribution of visible matter may be far removed from 
the physically motivated axes of the dark matter halo; a 
problem exacerbated by the fact that an LCDM dark mat- 
ter halo extends far beyond the baryons. One possible way to 
circumvent this degeneracy is to use a population of galaxies 
that have a known, preferred direction that can serve as the 
basis of a coordinate system that is more readily linked to 
the physical system of the galaxy and its dark matter halo. 

In this article we outline such a method, one which uses 
disk-like galaxies because the normal to the disk (which can 
be approximated using the projected galaxy shape) provides 
the necessary direction of preference, and apply it to data 
from SDSS DR6, the sixth data release of the Sloan Digital 
Sky Survey (lAdelman-McCarthv et al.l l2007h . We also use 
a parameterization of the angular distribution of satellite 
galaxies and present our results as an allowed region on an 
exclusion plot. This style of presentation extends previous 



studies which answer the question of whether or not satellite 
galaxies are distributed isotropically by addressing the ques- 
tion of what distributions are allowed. Our discussion will 
proceed with an outline of the formalism that we apply, our 
data selection criteria, our analysis of those data and its re- 
sults, and a discussion of the implications of our results, pos- 
sible ways to improve the sensitivity of this technique, and 
how this technique can complement other, existing methods 
used to observationally study dark matter halos. 

Before we continue, we note an apparent source of con- 
fusion in the literature regarding the coordinate systems 
used in numerical studies. In particular, there is a subtle but 
important difference between the moment of inertia tensor 
and the tensor of cartesian second mass moments. For ex- 
ample, the zz component of the moment of inertia tensor I 
is 



l„ = / p(x + y 2 )dV 



(i) 



while the zz component of the tensor of cartesian second 
mass moments M is 



M zz = / pz 2 dV. 



(2) 



where p is the mass density, and V is the volume occupied by 
the matter. We recommend that authors clearly state what 
quantities they use in their work and take care when using 
the relevant terminology since the largest moment of inertia 
corresponds to the smallest cartesian second moment and 
vice versa. In the remainder of this work we use the cartesian 
second moments (or just "second moments") exclusively. 



2 APPROACH AND FORMALISM 

Throughout this article we refer to the "physical" and the 
"projected" coordinate systems (see Figure [T}. The latter 
system is based upon the distribution of the luminous mat- 
ter where the primary axis is along the observed major axis 
of the galaxy. The former is the coordinate system based 
upon the angular momentum vector of the galaxy. Relat- 
ing this coordinate system to observation requires some as- 
sumptions, namely: 1) that the angular momentum vector 
of a galaxy is normal to the disk (or highly correlated with 
the normal) and 2) that the inclination of the galaxy can 
be found using the axis ratio of the isophotal contours. We 
will state additional assumptions as they occur. While not 
technically an assumption of our analysis, our language will 
imply that the angular momentum vector of the galaxy is 
correlated with a principle axis of the dark matter halo (pre- 
sumably either the largest for a prolate halo or the smallest 
for an oblate one). 

We proceed with the premise that if there is some 
anisotropy in the distribution of satellites, then the nor- 
mal to the disk is the natural axis about which to measure 
the deviations; and that the system can be approximated as 
cylindrically symmetric. That is, in the plane of the disk the 
effects of a triaxial dark matter distribution are negligible 
or are washed out due to randomness in the orientations of 
the hosts. Therefore, a sample of disklike galaxies of vari- 
ous projections allows us to constrain the allowed physical 
distributions of satellite galaxies. 

For a cylindrically symmetric (though highly 
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Figure 1. a) Diagram of the physical coordinate system con- 
structed about the angular momentum vector of the galaxy. The 
angle 8 is the polar angle in a spherical coordinate system, b) 
Diagram of the projected coordinate system constructed about 
the major axis of the observed luminous matter. The angle tp is 
the azimuthal angle in a cylindrical coordinate system where the 
axis of symmetry is along the line of sight. 

anisotropic) distribution of satellites, the projected 
distribution of satellites is a strong function of the inclina- 
tion of the galaxy with respect to the line of sight as shown 
in the top panel of Figure [2] The bottom panel of Figure [2] 
shows how the projected cumulative distribution function 
(CDF) differs from one line of sight to another. For the 
projected CDF, the angle tj) represents the minimum angle 
between the line connecting the centroids of the host and 
the satellite and the long axis of the projected disk of the 
host; the absolute value of ip must be between and n/2 
(see Figure [TJ). 

In order to study the physical CDF, we f irst utilize the 
param eterization of the distribution used by IZentner et all 
(2005), which is a power law of the form 

CDF(a;) = ax + (1 - a)x b (3) 

where x = | cos(0)| and where 9 is the polar angle of the 
physical coordinate system. In this parameterization, a char- 
acterizes the relative isotropy of the distribution of satellites 
and b encapsulates some functional form of the deviations 
from isotropy. Our analysis approach is designed to constrain 
the physical distribution of satellite galaxies using a set of 
projected CDF's, each of which corresponds to a different 
inclination angle of the galaxy with respect to the line of 
sight. 

By requiring that the physical CDF begin at zero, that 
it end at unity, and that it has a nonnegative derivative on 




tt/6 7i/3 n/2 

w 

Figure 2. Top: Projected distributions of a cylindrically symmet- 
ric and highly anisotropic, prolate distribution of satellite galaxies 
(grey dots) with different inclinations with respect to the line of 
sight. Bottom: The cumulative distribution function of satellites 
for each of the projections. The value of the angle ifi is the small- 
est angle between the long axis of the projected disk and the line 
connecting the centroids of the host and satellite and can have 
values between and w/2. If the distribution were oblate, the 
CDF would be consistently above the diagonal. 

the interval from to 1, the allowed values for the parame- 
ters a and b are restricted to lie in the unshaded region shown 
in Figure [3] While this parameterization has a straightfor- 
ward functional form, the infinite extent of the footprint -due 
to the hyperbolic boundary-makes it ill suited for an exclu- 
sion plot. Consequently, we transform the (a, b) parameter 
space to the (h, k) parameter space where the point (h, k) 
denotes the point where the physical CDF deviates maxi- 
mally from the isotropic distribution as shown in Figure [4] 
The transformations between the (a, b) representation and 
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Figure 3. Allowed values (unshaded) for the parameters a and b 
from equation [3] 



Physical CDF 




0.2 0.4 0.6 0.8 1 

|cos(0)| 

Figure 4. Illustration of the relationship between the parameters 
a and b from equation \3\ and the parameters h and k. The point 
[h, k) corresponds to the location where the physical CDF of the 
satellite galaxies deviates maximally from the isotropic case. 

the (h, k) representation are most concisely given by 

h = 6 1/(1 - 6) (4) 

and 



There is a degeneracy in this representation that corre- 
sponds to the isotropic case (where b — 1). In principle this 
may be an issue, but in practice any problem that arises can 
be avoided by choosing a representative point for that sce- 



Figure 5. Allowed values for the parameters h and k (unshaded) 
where the point (h, k) corresponds to the location of the point on 
the physical CDF which deviates maximally from the isotropic 
case. This representation allows only finite values for both pa- 
rameters. 

nario. Leaving that item aside, the allowed values of h and 
k are shown as the unshaded region in Figure [5] Our results 
will be shown as an allowed region on the (h, k) plane. 



3 ANALYSIS APPROACH 

In order to identify the values for h and k that are consistent 
with data, we divide the h — k plane into a 65 x 100 grid. 
For each point in the allowed parameter space (4290 points) 
and for each bin we did the following: 

(i) Generate a set of satellite galaxies derived from the 
point (h, k) whose physical distribution is uniform and ran- 
dom in azimuth and uniform but equally spaced in the phys- 
ical CDF-that is, for each, equally spaced value between 
and 1 we solve equation © for 8. This set is at least 15000 
times larger than the number of of satellite galaxies in each 
bin and serves as the population for our statistical analysis. 
The exact number for each bin (~ 7 x 10 B ) is chosen in order 
to avoid the need to interpolate between points during the 
analysis. 

(ii) Rotate the population by the midpoint angle of the 
corresponding bin and rewrite the resulting distribution in 
the projected coordinate system. 

(iii) From the projected population, draw 1000 samples of 
galaxies. The size of these samples are equal to the number 
of data in each bin. 

(iv) Calculate the sum of the squares of the differences 
(x 2 ) between the data CDF and the population CDF as well 
as the x 2 f° r each of the sample CDF's and the population 
CDF. 

(v) Count the number of samples whose \ 2 ls smaller than 
that of the data and assign this value to the point (h, k). We 
consider this value (divided by 1000) to be the probability 
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Figure 6. Allowed values for h and A; for the example shown in 
Figure [2] The contours correspond to the 90%, and 99% confi- 
dence levels. The large gray regions are where the values of h and 
k do not satisfy the requirements of a cumulative distribution 
function. 

that the data were drawn from the distribution associated 
with (h, k). 

After completing the above steps for each bin, we make 
a likelihood function by multiplying the three probabilities 
for each point (h, k). Finally, we identify the contours where 
the likelihood is greater than 0.1 and 0.01 (corresponding to 
the 90% and 99% confidence intervals). 

3.1 Example 

Prior to our analysis of SDSS data, we present an example 
in order to demonstrate the kind of result that this approach 
yields. Consider the distribution shown in Figure [2] The pa- 
rameters that we used to generate those plots are h = 0.8 
and k = 0.15 (a = 0.115, b = 12.2). We selected a sample of 
1200 satellites from that distribution-400 into each of three 
projection bins with boundaries of 0°, 36°, 54°, and 90°. 
Each satellite is selected as though it were viewed from a 
random orientation within its corresponding bin. The selec- 
tion properties and the binning of this sample are chosen to 
mimic those of the data as explained in the next section. 

We apply the analysis approach outlined above to this 
sample. The allowed values of the parameters h and k are 
shown in Figure |S] where the contours correspond to the 90%, 
and 99% confidence levels. The white dot is the location of 
the true values of the parameters that were used. 



4 HOST AND SATELLITE SAMPLE 
SELECTION 

We wish to apply the analysis approach from section [3] to 
a sample of real galaxies. The galaxies for this study, both 
hosts and satellites, are selected from the sixth data release 



of the Sl oan Digital Sky Survey spectro scopic and galaxy 
catalogs l|Adelman-McCarthv et al1l2007l ). We select galax- 
ies with redshifts between 0.005 and 0.3 and require that 
their redshifts have a confidence level (zconf ) greater than 
0.9. To select our sam ple of disk- l ike ga laxies, we use the 
colour divider given in lYang et all (|2006l ) 

01 ( fl -r)<0.83. (6) 

This criterion requires that we correct for extinction and 
that we make k corrections to a redshift of 0.1 (as indi- 
cated by the supers cript in equatio n ©); for which we use 
the kcorrect utility (Blan ton et al.ll2003T ). We include in the 
sample all galaxies where the axis ratio of the light distribu- 
tion (using the r-band isophotes) is less than 0.2 since these 
are likely to be disk galaxies reddened by dust. We also use 
the r-band isophotal axis ratio to determine the inclination 
of each galaxy with respect to the line of sight by assuming 
that 

cos(i) = B/A (7) 

where A and B are the angular extent of the major and 
minor axes respectively. 

We select our sample of candi date hosts galax ies and 
satellites following the criteria from iBrainerdl |2005i ). That 
is, a host galaxy must be brighter by one magnitude than 
any other galaxy: 1) within a projected circular apeture of 
radius 0.7Mpc and 2) whose relative velocity (calculated us- 
ing the spectroscopic redshift) is within lOOOkm/s of the 
host. The satellite galaxies are then selected in a similar 
manner. A satellite must be two r-band magnitudes dim- 
mer than its host, must lie within a circular apeture of ra- 
dius 0.5Mpc, and have a relative velocity within 500km/s. 
A host /satellites system is rejected if the total brightness of 
the satellites exceeds that of the host or if there are more 
than five satellites surrounding the host. These last criteria 
are to avoid cases where a single galaxy is deblended into 
multiple components. Our final sample includes 1279 host 
galaxies and 1595 satellites. 

We separate these galaxies into three bins of differing 
axis ratios such that the bins span the entire allowed range 
of inclinations between and 7r/2. We choose bins such that 
each of the bins have comparable numbers of satellites. The 
first bin is for host galaxies where the angle of inclination 
between the normal to the disk and the line of sight are 
between 0° and 36°-using equation (0, the second bin has 
inclinations between 36° and 54°, and the third bin has in- 
clinations between 54° and 90°. These bins contain 385, 514, 
and 380 host galaxies with 496, 628, and 471 satellite galax- 
ies respectively. Each bin has approximately 1.25 satellite 
galaxies per host. 

4.1 Analysis of SDSS Galaxies 

Figure [7] shows the projected CDF's for each of the three 
inclination bins. A Kolmogorov-Smirnov test between these 
different CDF's is inconclusive which indicates that they 
are consistent with being drawn from the same distribution 
(which implies that the satellites for this sample are nearly 
isotropic). By applying the procedure outlined in section [3] 
to these data we obtain the exclusion plot on the (h, k) plane 
that is shown in Figure [8] This figure shows that the most 
likely values for the parameters lie along the line h = k and 
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Projected CDF 
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Figure 7. The projected CDF's for the different inclination bins. 
The light gray points are for the highest inclination bin (nearly 
edge-on), the dark gray points are for the middle inclination bin, 
and the black points are for the lowest inclination bin. The line 
representing an isotropic distribution in the projected coordinate 
system is also shown. 
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h 

Figure 8. Allowed values of h and k for observed, blue galaxies 
from SDSS DR6. The contours correspond to the 90% and 99% 
confidence levels. 

indicates that the physical distribution of satellite galaxies 
surrounding disk-like hosts is nearly isotropic. The value of 
the liklihood function at its maximum is 0.80 (for the sample 
in section [4] it is 0.29) indicating that the model that we 
assume for the satellite galaxy distribution (|3]) is consistent 
with the SDSS data that we analyze. Moreover, it is unlikely 
that additional data could demonstrate the inconsistency of 
the model in the forseeable future. 



5 DISCUSSION 

Certain properties of our galaxy sample, our analysis ap- 
proach, and our results merit further discussion and as well 
as comparison with other methods employed to study the 
distribution of satellite galaxies. For example, one system- 
atic effect that we did not address is that disk galaxies are 
not infinitesimally flat , nor are the y perfectly axisymmetry 
when viewed face-on IjRvdenl 2006). The result is that the 
inclination of edge-on galaxies are underestimated by our 
approximation in equation while the inclinations of face- 
on galaxies are overestimated. 

The effect of these issues depends upon the physical dis- 
tribution of satellite galaxies, whether oblate or prolate, and 
upon whether the galaxy is edge-on or face-on. To explore 
the consequences of this systematic effect, consider four lim- 
iting cases: a maximally oblate distribution that is viewed 
edge-on, the same distribution viewed face-on, a maximally 
prolate distribution viewed edge-on, and a prolate distribu- 
tion viewed face-on. 

The maximally oblate case is where all the satellite 
galaxies lie in an annulus in the plane of the disk. If viewed 
edge-on, where the inclination is underestimated, the anal- 
ysis would assume that the apparent line of galaxies was ac- 
tually drawn from a distribution that was slightly rotated, 
which in projection would appear elliptical. Thus, the in- 
ferred physical distribution resulting from the analysis of 
edge-on galaxies is slightly more isotropic (less planar) than 
what is actually the case. For face-on galaxies with the same 
planar distribution, overestimating the inclination would in- 
correctly assume that the apparent cylindrically symmetric 
distribution of satellites was again viewed from a slightly 
different angle. In this case, the inferred distribution would 
have, incorrectly, more satellites along the major axis of the 
host. The effect of this on the overall inferred physical dis- 
tribution of satellite galaxies is less straightfoward than for 
the edge-on galaxies because of the assumption of cylindrical 
symmetry. Here the overestimated inclinations give an ob- 
served distribution of satellites that is always be anisotropic 
in the same manner; that is, about the major axis of the ob- 
served host. The result is an incorrect probability assigned 
to the parameters (h, k) for that projection bin-most of the 
parameters would have slightly lower probabilities than they 
should because of the violation of the cylindrical symmetry 
assumption. This effect is mitigated somewhat by the fact 
that the face-on galaxies contribute the least to the final 
liklihood function. 

For the maximally prolate case the satellites would lie in 
a narrow cone that opens above and below the host disk. For 
edge-on galaxies in this scenario, the inferred distribution of 
satellites would be less prolate (also more isotropic) than 
the actual distribution. For face-on galaxies the satellites, 
as with the oblate case for face-on galaxies, would no longer 
be cylindrically symmetric but would have more satellites 
near the observed minor axes of the host (this can be seen 
in the top panel of figure[2] which has different projections of 
a highly prolate distribution, the edge-on distribution would 
look less edge-on and the face-on distribution would look less 
face-on). Note that the effects for face-on galaxies would be 
relatively small and contained within the corresponding pro- 
jection bins. The more important point of this discussion is 
that for edge-on systems-which have the highest sensitiv- 
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ity to the physical distribution of satellites-the incorrectly 
assumed inclination results in an inferred satellite distribu- 
tion that is more isotropic than what is present in truth. 
This is true regardless of whether the distribution is oblate 
or prolate. 

Another issue to explore is the effect of binning the 
galaxies into either more or fewer bins. Previous studies ef- 
fecti vely used o ne bin , typically all galaxies with B/A < 0.8 
(e.g. IBrainerdl l|2005h ). While in principle it may be possi- 
ble to identify the physical distribution of satellites by this 
single, projected distribution (since they are both functions 
of one variable), there are some systematic effects that this 
study has brought to light which may complicate that task. 
For example, the host galaxies should be randomly oriented 
with respect to our line of sight. Yet, the data that we ana- 
lyze do not have this property. 

A likely cause for this discrepancy, particularly for 
highly inclined galaxies, is dust which affects bo th the color 
of a galaxy and its luminosity l)Shao et al.ll2007T ). Thus, the 
color cut that we apply does not select a uniform sample as 
a function of inclination since edge-on galaxies will appear 
more red. In addition, extinction causes the edge-on galax- 
ies to appear fainter than their face-on counterparts. Thus, 
one would expect fewer edge-on host galaxies in our sample. 
This inclination-dependent completeness would need to be 
addressed if one were using only a single bin of galaxies to 
infer the physical distribution of satellites. 

By increasing the number of inclination bins, one need 
not assume nor require that the hosts be uniformly sampled 
in inclination since each bin is largely independent of the 
others. Consequently, the results presented here suffer less 
from this effect than an attempt to identify the satellite dis- 
tribution from a single bin. There is no particular limit to 
the number of inclination bins that one might use to con- 
strain this distribution-we use three for largely computa- 
tional considerations. However, once the number of galaxies 
in each bin becomes too small (< 10) then it would likely be 
better to alter the approach and use the PDF rather than 
the CDF, identifying the probability that a particular satel- 
lite would be found near a particular angle given the host 
inclination and the distribution that is then being tested. 
Another consideration is that if the difference in inclination 
between bins is less than the uncertainty in the estimates of 
the inclination, then there would be a significant amount of 
overlap between the members of different bins which would 
complicate the analysis and possibly the interpretation of 
the results. 



6 CONCLUSIONS 

Our results have implications for the distribution of matter 
in dark matter halos, the orientation of the disk inside of 
a dark matter halo, and (or) dynamical effects that depend 
upon the mass of the satellite galaxy. Whi l e thes e results 
are consistent with the results of Yang et al.l (|2006h for blue 
host galaxies, there may be some reasons that observed dis- 
tribution might erroneously imply an isotropic one. These 
include: 1) the axis ratios of the disk-hosting dark matter 
halos might be near unity, resulting in a nearly spherical dis- 
tribution of satellite galaxies such that our sample cannot 
distinguish between the actual distribution and a spherical 



one-something that may be resolved with a larger sample of 
satellites; 2) the galaxy disks may not be highly correlated 
with the major or minor axes of the dark matter halo-if 
the normal to the galactic disk is commonly oriented > 20° 
away from the axes of the dark matter halo then the dif- 
ferent projection bins would mix and any signal from the 
dark matter halo would be blurred ( see iBailin fc Steinmet j 
l|2005l ); [Agust sson fc Brainerdl i|2006n ); 3) the selection cri- 
teria that we use may not provide an unbiased sample of 
all of the satellite galaxies in a system-perhaps the dynam- 
ics of the relatively bright satellites that we use (ones that 
would be selected for the SDSS spectroscopic survey) trace 
the dark matter halo less strongly than smaller galaxies that 
were not selected. 

Future applications of this technique may be able to ad- 
dress some of these issues. Given the fact that the SDSS data 
contain millions of galaxies and we only analyze ~ 1600, 
there may be a more sophisticated selection criteria that 
would identify a larger number of satellite galaxies, espe- 
cially smaller satellites whose dynamics are more strongly 
dominated by the host. In particular, if a method were 
found that could use photometric redshifts instead of spec- 
troscopic redshifts in the selection criteria, then a much 
larger sample could be drawn; this would be very useful 
for data from other existing and planned photometric sur - 
veys (e.g. Dar k Energy Su rvey (|DES C ollaboration 2007), 
PanSTARRS (|Kaiserll2002h . or LSST (|LSST Collaboration! 
12007ft ). In addition, other selection criteria may eliminate or 
red uce systematic effec ts that could cause a biased sample 
fsee lBailin et alJ (|2007h ). We chose our selection criteria as a 
means to tie this technique to previous studies. Implement- 
ing a different or more sophisticated selection criteria lies 
outside the scope of this work. 

If none of the above issues aff ect our s a mple then 
the observed ani sotropy reported in IBrainerdl (2005) and 
lYang et al.l l|2006h could be explained as the sum of a spher- 
ically symmetric distribution around blue galaxies and an 
anisotropic distribution around red gal axies. This explana - 
tion is consistent with the findings of lYang et al.l (2006). 
The anisotropic distribution of satellites around red galax- 
ies may imply that the merging history plays a major role 
(if not the primary role) in the triaxiality of dark matter ha- 
los and the resulting distribution of satellite galaxies should 
they prove an unbiased tracer of the mass when the sample 
is sufficiently large. 

Another way to confirm the conclusion that disk-like 
galaxies reside in nearly spherical halos is by applying this 
technique to gravitationally lensed systems-both weakly and 
strongly-where the host galaxy acts as the lens. While chal- 
lenging to implement, these probes characterize all of the 
matter in the lensing system and, in conjunction with the 
image of the lensing galaxy itself, could help identify the rel- 
ative orientation of the baryons and the dark matter halo. 
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